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ABSTRACT 

We outline a new method to derive a "snapshot" metaUicity estimate of stellar systems (providing 
one resolves at least the brightest part of the CMD) just on the basis of low-resolution (i.e., 6- 
8 A FWHM) spectroscopy of a small stellar sample. Our method relies on the Fe5270 Lick index 
measurements and takes advantage of the special behavior of this spectral feature, that reaches its 
maximum strength among the ubiquitous component of K-type giants. This makes the Fe5270niax 
estimate a robust and model- independent tracer of cluster [Fe/H], being particularly insensitive to the 
age of the stellar population. 

A comparison of the Fe5270max distribution derived from globular and open clusters, as well as from 
the field giant population in the Galaxy disk, confirms a tight correlation of the index maximum vs. 
cluster [Fe/H] allover the entire metaUicity range for stellar population with [Fe/H] > —2.0. Relying 
on a theoretical calibration of the feature, we trust to effectively infer cluster metaUicity within a 
typical uncertainty of 0.1-0.2 dex, depending on RGB luminosity sampling of the observations. 

A handful of stars (5-10 objects) is required for the method to be applied, with low-metallicity 
stellar populations more easily managed, being Fe5270max located within the few brightest RGB stars 
of the system. In any case, we show that even the observation of a coarse stellar set would allow us 
to place a confident lower limit to cluster metaUicity. 

Subject headings: stars: abundances — globular clusters: general — open clusters and associations: 
general — Galaxy: disk 



1. INTRODUCTION 

A fair estimate of metal abundance in stellar systems 
still remains a central issue for any detailed assessment of 
the other cluster distinctive properties. MetaUicity mod- 
ulates, in fact, both effective temperature and apparent 
colors of stars, while the internal composition affects the 
nuclear engine, and therefore stellar lifetime. Even in 
case of resolved stellar systems, this entangled behavior 
may lead to a biased interpretation of cluster age, based 
on the CMD morphology, an effect often referred to as 
the "age-metallicity dilemma" (Renzini & Buzzoni 1986; 
Worthey 1994). Such induced age uncertainty also re- 
flects in the distance determination of star clusters, as 
far as one tries to compare the apparent magnitude of 
the Main Sequence Turn Off (TO) point with the ap- 
propriate theoretical luminosity to derive therefrom the 
distance modulus. One further difhculty also deals with 
the proper assessment of dust reddening, that may affect 
CMD morphology leading, in general, to an artificially 
enhanced value of [Fe/H].-'^ 

To overcome these problems, one would like to prelim- 
inarily derive the cluster metaUicity from accurate abun- 
dance analysis of individual stars (typically red giants, 
due to their brighter intrinsic luminosity at optical wave- 
lengths) through high-resolution spectroscopy (e.g., Car- 
retta & Gratton 1997; Kraft & Ivans 2003). This dehcate 
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^ In terms of color excess, this effect can be quantified in E{B — 
V) ~ 0.14 A [Fe/H] (Buzzoni 1995a). 



task, however, is extremely time consuming, and a far 
more straight shortcut is often pursued relying on inte- 
grated cluster spectroscopy, usually taken at much lower 
resolution. Narrow-band spectrophotometric indices, as 
derived from integrated low-res (~ 6-8 A FWHM) ob- 
servations, usually provide the basic diagnostic scheme 
to match theoretical models from stellar population syn- 
thesis and derive therefrom cluster properties (e.g., de 
Freitas Pacheco et al. 1998; Beasley et al. 2002; Strader 
& Brodie 2004). The Lick system (Burstein et al. 1984; 
Worthey et al. 1994; Trager et al. 1998) stands out as the 
widest and most popular reference, assuring a systematic 
coverage of the main spectral features across the 4100- 
6400 A wavelength range, easily accessible from ground- 
based observations. 

Besides the advantage of this strategy, one can still 
question whether it effectively allows us to disentangle 
any age-metallicity degeneracy. In fact, integrated in- 
dices are sensitive both to intrinsic elemental abundance 
and to the temperature distribution of the underlying 
stellar population, thus delivering a composite and likely 
non-univocal piece of information (see a discussion in 
Buzzoni 1995b; Worthey et al. 1995; Tantalo & Chiosi 
2004). 

An alternative, and possibly more proficient way out, 
at least for resolved stellar clusters, can be envisaged 
and will be explored in this paper. It relies on a minimal 
"tuned" sampling of individual red giant stars of a clus- 
ter, observed at low spectral resolution such as to derive 
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the maximum strength for the Fe5270 index, which is by 
far the most popular one among the several Fe I features 
included in the Lick system. As we will see in Sect. 2, 
for their special physical properties, Fe i features prove 
to be ideal candidates for a reliable "snapshot" estimate 
of stellar system metallicity. 

Our method complements previous studies in the same 
line. For example, by relying on low-resolution near- 
IR observations of the strong Ca ii triplet among RGB 
stars, Da Costa & Hatzidimitriou (1998) and Parisi et al. 
(2009) recently recovered the [Fe/H] value for the SMC 
and its cluster system within less than a 0.2 dex uncer- 
tainty compared, for instance, to the high-res analysis of 
Battaglia et al. (2008). 

In the following discussion, we will outline our method 
in Sect. 3, and propose a plain theoretical calibration 
based on population synthesis models. Theoretical pre- 
dictions will be compared with a grid of Galactic globular 
and open clusters spanning the whole range of metallic- 
ity. In Sect. 4 we will summarize our conclusions and also 
briefly outline further applications of the method, also in 
view of the deeper forthcoming surveys of external stel- 
lar systems (i.e., extragalactic star clusters around and 
within nearby galaxies in the Local Group, etc.) taken 
with the new-generation telescopes. 

2. NARROW-BAND INDICES OF STARS 

Extended grids of synthetic stellar spectra, from which 
to derive Lick-index behavior vs. logTeff, logg, and 
[Fe/H], have been computed, among others, by Mould 
(1978), Barbuy (1994), Tripicco & Bell (1995), Chavez 
et al. (1996), Barbuy et al. (2003), and Bertone (2002). 
In addition, semi-empirical calibrations based on more 
or less complete samples of disk stars are due, among 
others, to Buzzoni et al. (1992), Buzzoni et al. (1994), 
Gorgas et al. (1993), Worthey et al. (1994), Borges et al. 
(1995) and Franchini et al. (2004). The latter results are 
often summarized in analytical sets of so-called "fitting 
functions" . 

Both approaches resulted in a complementary contri- 
bution as far as the physical behavior of the different 
spectral features vs. stellar fundamental parameters is 
concerned. In this respect, a recognized advantage of 
synthetic spectra is to directly explore index changes 
within a "controlled" grid of input parameters. However, 
models are still largely inadequate in their physical treat- 
ment of cool stars, (i.e., Tes ^ 4000 K, including late-K 
and M spectral types; see Bertone et al. 2004, 2008), 
when convection and the increasing effect of molecular 
opacity severely perturb the atmosphere structure and 
the emerging spectral energy distribution of stars at op- 
tical wavelength. 

From its side, the empirical fitting-function technique 
has the advantage of giving a simple and accurate sum- 
mary of index behavior across the full parameter domain 
of real stars. This partly overcomes any modeling uncer- 
tainty, and the analytical set of equations can easily be 
implemented in theoretical models of stellar populations 
to compute the integrated indices of a system without re- 
quiring any spectral input for individual stars (Buzzoni 
et al. 1992, 1994; Worthey 1994; Maraston et al. 2001; 
Maraston 2005). As a drawback of this approach, how- 
ever, one has to remark the evident limitation of any em- 
pirical star sample (forcedly dominated by a component 
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Fig. 1. — The Fe5270 index as a function of T^g from Bluered 
synthetic stellar spectra of model atmospheres with log </ = 3.0 dex 
and solar metallicity (dots). Square markers trace the raw Fe5270 
index trend by neglecting the TiO molecular opacity. 

of solar neighbors) to properly probe the full metallicity 
scale outside the solar range. 

2.1. Fe5270 diagnostic properties 

Among the set of Lick indices in the optical wavelength 
range, the Fe i lines reveal potentially the most promising 
ones for the present analysis. This is especially true for 
the Fe5270 feature, that begins to appear in stars cooler 
than Tcs — 7500 K and reaches a maximum equivalent 
width about 4200-4600 K (Buzzoni et al. 1994). The lat- 
ter is the typical temperature of faint low-main sequence 
dwarfs and, most importantly, also of bright MK III gi- 
ants, which populate the AGB/RGB region of the CMD 
diagram. 

Although not so prominent, the Fe5270 feature (and 
its close neighbor Fe5335) can easily be recognized also 
in the integrated spectra of elliptical galaxies (Worthey 
et al. 1992; Trager et al. 1998), and its presence is known 
to trace the relative contribution of low- RGB stars in the 
galaxy stellar population Buzzoni (1995a). 

In Fig. 1 we explored in finer detail the Fe5270 diag- 
nostic properties from the theoretical point of view, by 
means of an ad hoc grid of Bluered synthetic spectra 
(Bertone 2002; Bertone et al. 2008). Spectral synthe- 
sis relied on Synthe code (Kurucz 1993; Sbordone et 
al. 2004) and the grid of Atlas9 revised model atmo- 
spheres (Castelli & Kurucz 2003) along the 3500-4500 K 
Tcs range. For Tes ^ 4000 K, the main contributor to 
line blanketing, along the three wavelength bands that 
define the index, is found to be Fe i, with a fine structure 
of more than 200 lines; particularly strong is, of course, 
the A5269.5 line. However, in this temperature regime, 
also important is the contribution of neutral Titanium 
and Chromium and, secondly. Calcium and Cobalt. 

The TiO opacity becomes increasingly stronger with 
decreasing temperature, and becomes dominant for 
Toff < 3750 K. As a result, the prevailing effect of the 
TiO molecular bands on the Fe5270 index is in the sense 
of depressing the red-side pseudo-continuum, thus reduc- 
ing the apparent strength of the index. According to the 
models, one may even conclude that the TiO molecule 
is actually the main responsible for constraining in tem- 
perature the Fe5270 peak value (see, again. Fig. 1). 

We can probe index sensitivity also by means of empir- 
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Fig. 2. — Fitting-function predictions for Fe5270max vs. [Fe/H] 
relationship according to Buzzoni et al. (1994) (B94, see eq.l), 
Worthey et al. (1994) (W94, eq.3) and Gorgas et al. (1993) (G93, 
eq.4). Two relevant cases for stellar gravity are considered, namely 
log g = 3 and 1 dex. The theoretical values from Bluered synthetic 
spectra are also overplotted (dashed line). 

ical data. From the Buzzoni et al. (1994) fitting function 
we can write in fact: 

Fe5270max = 1-15 [Fe/H] + 0.251og.g + 3.80, (1) 

that fairly well compares with the Bluered theoretical 
prediction,'^ namely 

Fe5270n,ax = (1-68 - 0.1 logg) [Fe/H] + 4.80. (2) 

According to Worthey et al. (1994) models we have 

Fe5270max = 1-29 [Fe/H] + (0.11 logg - 0.51) log5 + 4.54, 

(3) 

while Gorgas et al. (1993) results provide'^ 

Fe5270max = (1-79 + 0.05 logg + 0.19 [Fe/H]) [Fe/H] + 
+ (0.13 logg- 0.31) log5 + 4.55. (4) 

Note from the equations the weak index sensitivity to 
stellar surface gravity. For example, any change of one 
dex in logg only contributes at most by ±0.3 A to the 
Fe5270niax variation. In terms of overall properties of a 
stellar aggregate, this figure reflects in a change of TO 
stellar mass of roughly 0.3 dex, which implies a vari- 
ation of nearly one order of magnitude in stellar life- 
time. Therefore, as a major conclusion, we have that 
Fe5270niax in a stellar population is nearly independent 
from age. 

A comparison of the different calibrations is shown in 
Fig. 2, where we report the output index for two refer- 
ence values of stellar gravity (namely, log 5 = 1 and 3 
dex) appropriate for red giants. Within a remarkably 

^ The [Fe/H] dependence of the Bluered synthetic index has 
perhaps to be considered as an upper limit, as molecular absorption 
in the relevant wavelength range should likely be stronger than 
predicted by models at solar and super-solar metallicity. According 
to our experiments, this makes the current Bluered output to 
overestimate by ~0.5 A the real Fe5270 index strength as observed 
for metal-rich stars. 

3 Note that the Worthey et al. (1994) and Gorgas et al. (1993) 
results are not completely independent as both rely on a largely 
overlapping stellar sample. One relevant difference, however, is 
that Gorgas et al. (1993) use the V — K color instead of T^g for 
their fitting functions. 
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Fig. 3. — The expected Fe5270 Lick index strength in SSPs 
for an illustrative set of theoretical CMDs, from Buzzoni (1989, 
1995a) models. Left panel explores the index change vs. [Fe/H] 
(as labeled on the plots) for a set of 15 Gyr populations, while 
right panel takes into account the index variation for fixed values 
of metallicity with changing SSP age (from 15 and 2 Gyr). The 
width of the isochrones is proportional to the local index value. 
Small square markers locate the maximum value, throughout. 

good agreement among the different datasets, compared 
to Buzzoni et al. (1994) calibration, one notes from the 
plot that Worthey et al. (1994) fit displays an even lower 
dependence on log 5, while Gorgas et al. (1993) fit stands 
out for its steeper trend with metallicity, especially at 
super-solar regimes. Overall, we can conclude that the 
direct Fe5270 dependence on Fe abundance turns out to 
be AFe5270 = a A [Fe/H], with a = 1.2 ^ 1.4. 

3. TRACING STELLAR POPULATIONS 

In addition to its direct sensitivity to [Fe/H], a ma- 
jor supplementary advantage of using Fe5270 feature as 
metallicity tracer in stellar populations resides in the spe- 
cial property of the index to always peak within the tem- 
perature range sampled by red giants in stellar systems 
richer than [Fe/H] > —2 dex.** 

Conversely, this is not the case, for instance, of other 
stronger features, used to derive popular indices like the 
Lick Mg2 or the near-IR CaT triplet index (Jones et al. 
1984; Idiart et al. 1997; Cenarro et al. 2001). Both these 
features peak, in fact, among giants and dwarfs at the 
coolest temperature tail of each CMDs, and their depen- 
dence on SSP metallicity is, therefore, naturally prone to 
more entangled age effects and to other model-dependent 
features, such as convection and mass loss properties 
(both constraining the AGB/RGB tip location). 

The Fe5270niax location for the illustrative case of dif- 
ferent simple stellar populations (SSPs), with varying 
either age or [Fe/H], is shown in Fig. 3. A glance to 
the figure makes clear that Fe5270 homogeneously probes 
the ubiquitous component of K1-K3 HI stars within each 
stellar system. 

The collection of empirical estimates of Fe5270inax for 
several open and globular clusters of the Galaxy, accord- 
ing to the Gorgas et al. (1993) database, is summarized 
in Fig. 4. These results are collected in Table 1, where 
we report for each cluster the strongest value reached by 
the Fe5270 index among the sampled stars, together with 

* Apart from Fe5270, two other Lick indices display a similar 
behavior: the G-band and the CN molecular indices (see Gorgas et 
al. 1993). For its atomic origin, however, Fe5270 should in principle 
be a more accurate and unbiased metallicity tracer. 
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TABLE 1 

Summary of relevant data for Galactic open and 
globular clusters^"' 



Name 


TypeC') 


[Fc/H] 


Fe5270max 
fAl 


No. of 

sa,rriplGcl stars 


NGC 188 


O 


-0.02 


4.03 


24 


NGC 7789 


O 


-0.08 


4.45 


12 


M67 


o 


-0.10 


3.89 


20 


M7l('=' 


G 


-0.73 


3.18 


36 


M5 


G 


-1.27 


2.94 


5 


MIO 


G 


-1.52 


2.32 


4 


M13 


G 


-1.54 


2.09 


3 


M3 


G 


-1.57 


2.57 


3 


M92 


G 


-2.28 


1.43 
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Trom Gorgas et al. (1993) Table 3. Metallicity of globular clus- 
ters is from Harris(1996) , while for open clusters it derives from 
the WEBDA database (Paunzen 2008); 

— Open cluster; G — Globular cluster; 

'^A supplementary entry exists for this cluster with Fe5270max — 
3.99 but star has been dropped by Gorgas et al. (1993) as possible 
field interloper. 



the sample size in each system, and the adopted cluster 
metallicity. Note that, as Fe5270 strength peaks among 
increasingly brighter red giants with decreasing metallic- 
ity (see left panel of Fig. 3), just a handful set of bright 
stars in metal-poor globular clusters is sufficient to suit- 
ably pick up Fe5270max- A slightly deeper spectroscopy 
is required, instead, for metal-rich systems, but even in 
this case a dozen of stars taken at low spectral resolution 
are fully sufficient to suitably constrain the index tip. As 
an interesting example, in this regard, we also overplot 
in Fig. 4 the full sample of field giants (class MK III or 
1 < log 5 < 3) studied by Buzzoni et al. (1994), and Buz- 
zoni et al. (2001); it is evident that the upper envelope 
of the Fe5270 distribution for the entire stellar dataset 
effectively matches the overall cluster trend. 

The data of Fig. 4 are also compared with theoreti- 
cal predictions from population synthesis models. The 
displayed models refer to 12.5 Gyr SSPs with different 
metallicity from Buzzoni (1989, 1995a). The case of 
younger (5 Gyr) populations is also accounted for the 
metallicity range of open clusters.^ 

4. RESULTS AND CONCLUSIONS 

A fit to all synthesis models of Fig. 4, excluding the 
[Fe/H] — —2.27 case, provides 



Fe5270n 



= 1.40 



(±0.03) 



[Fe/H] + 4.52 



(±0.02)i 



(5) 

with the displayed coefficient uncertainty at la level. 
Compared to the original stellar fitting function of 
eq. (1), we report here a slightly steeper slope (i.e., a 
more sensitive dependence) in the Fe5270 vs. [Fe/H] SSP 
relationship (namely a = 1.40 here vs. 1.15 of eq. 1). 
This is due to the reinforcing effect of the log g variation, 
which increases among metal-rich giants, for fixed SSP 
age. 

^ For young (t < 6 Gyr) and extremely metal-poor ([Fe/H] < —2) 
SSPs, however, models indicate that red giants are "too blue" (i.e., 
too warm) for the Fe5270 feature to span its full range and reach 
its maximum. In these case the Fe5270 index mat actually drop 
to nominal values and even remains undetected in the spectra of 
these stars. 
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Fig. 4. — The observed distribution of Fe5270max from the star 
cluster sample of Gorgas et al. (1993). Big dots mark the strongest 
value of Fe5270 among the sampled giant stars in each open and 
globular cluster, as summarized in Table 1. Only bona fide cluster 
members are considered, according to Gorgas et al. (1993) classi- 
fication. The two field giants samples studied by Buzzoni et al. 
(1994, 2001) (small triangles and squares, respectively) are also 
overplotted, considering only giant stars (i.e., 1 < logg < 3) with 
confident fundamental parameters. Observations are compared 
with stellar fitting functions for fixed gravity, i.e. logg = 3 and 1 as 
labeled, according to eq.(l). Also theoretical predictions from stel- 
lar population synthesis are reported, according to Buzzoni (1989, 
1995a) SSP models (thick vertical bars); an age t = 12.5 Gyr is 
adopted throughout, with the additional cases of 5 Gyr popula- 
tions for the metal-rich ([Fe/H] > —0.5) regime more pertinent to 
Galactic open clusters (left "column" at each relevant metallicity). 

Entering eq. (5) with the cluster values for Fe5270inax 
from Table 1 we could assess the statistical performance 
of our theoretical calibration as [Fe/H] predictor. One 
has to bear in mind, however, that this procedure is 
clearly biased toward underestimating cluster metallicity 
as we might be missing the real Fe5270 maximum. Of 
better interest may rather be to study the distribution of 
positive and negative residuals (in the sense "observed - 
predicted" [Fe/H]), separately. In summary, the relevant 
figures are the following: 



'^[Fe/H] = ±0-21 CICX 



^[l^e/H] = 0-13 dex 



^[Fe/H] = 0-26 dex 



(6) 



where the rms for the total sample residuals (l.h. side of 
the equation) comes in fact from a skewed distribution 
of positive (cpj^yjj]) and negative {c^^^yi]) residuals. For 
the previous arguments, the o'[fo/h] "^^-lue may actually 
provide a more genuine estimate of the real method per- 
formance, once a fair spectroscopic sampling of cluster 
stellar population can be assured. 

In conclusion, the claimed advantage of our diagnos- 
tic technique for a "snapshot" estimate of metal abun- 
dance in stellar systems (providing one resolves at least 
the brightest part of the CIVID) can be briefly summa- 
rized as follows: 

1) As Fe5270niax is negligibly affected by stellar grav- 
ity, its value within a stellar population is virtually age 
independent, and may therefore provide an unbiased 
measure of metallicity within a 0.1-0.2 dex internal un- 
certainty. The Fe5270max vs. [Fe/H] correlation is ex- 
pected to strictly hold for every stellar population with 
[Fc/H] > -2.0. ' 

2) As far as we know, within this framework, the 
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Fc5270 index performance is definitely better than that 
of any other stronger feature. For instance, by rely- 
ing on Mg2 or CaT* indices (see Cenarro et al. 2001), 
our procedure delivers cluster metallicity for the sam- 
ple of Table 1 with a minimum^ internal accuracy of 
c7([Fe/H])Mg2 0.25 dex and cr([Fe/H])caT* 0.42 dex, 
respectively. 

3) Compared to other spectroscopic methods, ained 
at deriving cluster metallicity from low-resolution spec- 
tra of their bright stars, our method has the superior 
advantage that it does not require any a priori assump- 
tion/inference about stellar temperature and gravity. It 
is distance and reddening independent, and uses an Iron- 
dominated feature (instead of any a-elcment feature, as 
for Mg and Ca indices) to infer [Fe/H]. 

4) Low-resolution spectroscopy {R = 600 800) is 
required for the method, planning just a "one-shot" ob- 
servation of the brightest stars of a given cluster with a 
standard multi-object (MOS) instrumental setup. Com- 
pared to broad-band imaging, one could envisage to de- 
rive useful spectroscopic information for stars about 3-4 
mag brighter than the imaging limit. 

5) As an added value, in addition to the Fe5270 index, 
spectra would naturally deliver supplementary informa- 
tion on other popular Lick indices, like Mg2, H/3, etc., at 
close wavelength. 

® Because one should also account for an additional uncertainty 
in case of unknown cluster age) 



6) A minimal set of giant stars is required (5 10 ob- 
jects), hopefully spanning the whole luminosity range. 
Low-metallicity stellar populations would be the favorite 
ones, as Fe5270max could easily be located among the 
very few outstanding stars of the cluster. Deeper obser- 
vations (down to the base of the RGB) would be required, 
instead, for metal-rich populations. According to Ren- 
zini & Buzzoni (1986), however, the luminosity-specific 
full number of Post-MS stars in a SSP with TO mass, 
Mto, is 

y~0.03MTo^^- (7) 

This means, for instance, that by sampling some 10^ Lq 
in the uncrowded external region of a globular cluster one 
may successfully pick up Fe5270niax among the few dozen 
brightest stars in the field. In any case, the observation 
of just the few outstanding members would confidently 
place a lower limit to cluster metallicity in lack of any 
other piece of information. 

It is a pleasure to thank Guy Worthey, the referee of 
this paper, for his valuable suggestions and for an ac- 
curate review of our work. We also acknowledge partial 
financial support from Mexican CONACyT under grant 
49231-E. 
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